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We present a fundamental plane (FP) analysis of 141 early-type galaxies in the Shap- 
ley supercluster at z=0.049 based on spectroscopy from the AAOmega spectrograph 
at the AAT and photometry from the WFI on the ESO/MPI 2.2m telescope. The 
I key feature of the survey is its coverage of low-mass galaxies down to (T~50kms~^. 

We obtain a best-fitting FP relation reOccr^ °6±° °^(/);:° *^^^° °^ in the R band. The 
shallow exponent of do is a result of the extension of our sample to low velocity dis- 
persions. Considering the subsample of do > 100 km s~^ galaxies, the FP relation is 
pr^ ' reOC(7^ '^^(/)g ^ '^^j consistent with previous studies in the high-luminosity regime. We 

OO ! investigate the origin of the intrinsic FP scatter, using estimates of age, metallicity and 

CO ■ a /Ye. We find that the FP residuals anti-correlate (>3cr) with the mean stellar age 

' in agreement with previous work. However, a stronger (>4(t) correlation with a/Fe is 

, also found. These correlations indicate that galaxies with effective radii smaller than 

' those predicted by the FP have stellar populations systematically older and with a 

\ over-abundances larger than average, for their a. Including a/Fe as a fourth parameter 

. in the FP, the total scatter decreases from 0.088 dex to 0.075 dex and the estimated 

intrinsic scatter decreases from 0.068 dex to 0.049 dex. Thus, variations in a/Fe ac- 
count for ~30% of the total variance around the FP, and ~50% of the estimated 
■ intrinsic variance. This result indicates that the distribution of galaxies around the 

?-H ' FP are tightly related to the enrichment, and hence to the timescale of star-formation. 

Our results appear to be consistent with the merger hypothesis for the formation of 
ellipticals which predicts that a significant fraction of the scatter is due to variations 
in the importance of dissipation in forming merger remnants of a given mass. 

Key words: galaxies: abundances, galaxies: ellipticals, galaxies: formation, galaxies: 
fundamental parameters, galaxies: structure 



1 INTRODUCTION 

Early-type galaxies are observed to obey a set of scal- 
ing relations that connect their photometric and kine- 
matic properties (e.g. Kormendy relation, Kormendy 1977, 
Faber- Jackson relation, Faber & Jackson 1976). Among 
these, the most notable, due to its surprising small scat- 
ter (~0.1dex), is the relation between effective radius r^, 
mean surface brightness within the effective radius (J)^ 



* E-mail: gargiulo@na.astro.it 



and central velocity dispersion ctq (|Diorgovski fc Davislll987l : 
iDressler et al.|[l987l ). In the three-dimensional space (logre, 
log (To, log(7)j,), elliptical galaxies populate a tight plane 
known as the fundamental plane (FP) and usually expressed 
in the form: 

logre = alogcro -|-/31og{/)p -f 7. (1) 

If elliptical galaxies formed a homologous family, i.e. systems 
with density, luminosity and kinematical structures equal 
over the entire early-type sequence and with constant mass- 
to-light ratios, then the virial theorem predicts a correlation 
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with a=2, 13=—1. However, observations show that the plane 
is somewhat "tihed" with respect to viri al expectations, 
with best-fit scalings a~1.3, 0.8 (e.g. Ij0rgensen et ahl 
1 19961 ). 

The origin of the FP tih has been much debated and 
can be interpreted as the breakdown of either of the two 
assumptions in the virial expectation. A systematic vari- 
ation in the mass-to-light ratio along the FP could be 
due to variations in the stellar content (age, metallicity 
or IMF) and/or the amount of dark matter among ellipti- 
cals. Performing detailed dynamical analyses of 25 galax- 
ies with SAURON in tegral-field stellar kinematics to r^, 
ICappellari et al.1 (|2006l ) find the "tilt" almost exclusively due 
to real M/L variations of the form (M/L)ocM°-^^='=°"^, while 
structural and dynamical non-homologies have negligible ef- 
fects. They also find the variation of the dynamical M/L ra- 
tio to correlate with the H/3 line-strength, and ascribe most 
of the tilt to stellar popula tion (a ge) effects. On the other 
hand (e.g. La Barbera et al. l2008al .Truiillo et al. l2004 l other 
authors find that the tilt is not primarily driven by stellar 
populations, but instead results from other effects, such as 
non-homology. 

Although the FP relation is quite tight, there is none 
the less a significant scatter around the plane that cannot 
be attributed to measurement errors. The origin of this in- 
trinsic component has been investigated by many authors. 
|j0rgensen et al.l (|l996l l found that they were unable to re- 
duce the scatter by introducing additional parameters, such 
as ellipticity or isophotal shape of the galaxies, into the FP 
relation. Variations in stellar populations along the sequence 
of early-type galaxies are found to be partially responsible 
for the intrinsic scatter (e.g. Gregg 1992; Guzman & Lucey 
1993; Guzman, Lucey, Bower 1993). Prugniel and Simien 
(1996) studying the correlation between the residuals from 
the FP and the residuals from the colour and Mg2 line- 
strength vs. luminosity relations, found that blue and low- 
Mg2 elliptical galaxies deviate systematically from the value 
predicted by the FP. Following this evidence Forbes et al. 
(1998), studying a sample of non-cluster galaxies, found that 
the residuals of the FP correlate with the ages of the galax- 
ies, i.e. that the scatter of the FP is partly due to variation 
in galaxy age at a given mass, and in particular to vari- 
ations in the time of the last starburst. On the contrary, 
they found that the effect of changes in metallicity is neg- 
ligible. Similar results were obtained by Reda et al. (2005) 
analysing a sample of isolated galaxies: some objects devi- 
ate from the FP relation having lower M/L ratio and this 
was interpreted as due to their younger stellar populations, 
probably induced by recent gaseous merger. The same con- 
clusions were reached by Wuyts et al. (2004) for two high- 
redshift clusters. They found that the residuals from the FP 
correlate with the residuals from the H/3 - ao relation. This 
confirms the role played by stellar populations in determin- 
ing the appearance of the FP, with relations appearing more 
dispersed for samples of galaxies that are more dispersed in 
age. 

The existence of the FP, its small observed scatter and 
the tilt have presented a long standing challenge to theo- 
retical models explaining the origin of early-type galaxies. 
In fact, whatever the scenario of formation and evolution of 
early-type galaxies is, it has to be able to explain the exis- 
tence of such a tight correlation and its deviation from virial 



expectations and therefore to link galaxy structure and dy- 
namics with their star-formation histories. 

In the recent years, observations and simulations have 
broadly supported the galaxy merging scenario which fits 
naturally into the ACDM hierarchical cosmology (e.g. 
Steinmetz & Navarro, 2002). In the hierarchical sce- 
nar io ellipticals form through the mergin g of disk galax- 
ies (|Toomre fc Toomre|[l973 : lToomrelll977^ . In the merging 
context, significant new insights have been made through 
large-scale gas dynami cal simulations of galaxy mergers 
iRobertson et all 120061 ). indicating that for lower mass 
galaxies, dissipation becomes increasingly important, driv- 
ing nuclear starbursts that contribute larger mass fractions, 
and producing systematic trends with mass in both the 
structures and stellar populations of the remnant ellipticals 
iHopkins et al.ll2008l ). The FP tih then arises as a direct 
consequence of the systematic trends with mass of the im- 
portance of dissipation during mergers. In the same scenario 
of galaxy formation, the origin of the intrinsic scatter in 
the FP arises as a combination of the scatter in the total 
baryon-to-dark matter content of the progenitor galaxies, 
and variations in the dissipational fractions at fixed stellar 
mass. This latter factor should be observable as correlations 
between the residuals from the FP and the stellar popula- 
tion para meters, and represen ts a critical test of the merger 
scenario (|Hopk ins et al.i2008l ) , through the predicted coevo- 
lution of the stellar populations and structures of elliptical 
galaxies. 

The increasing importance of dissipation in the forma- 
tion of low-mass galaxies and the different mechanisms that 
drive the evolution and the star-formation histories for low- 
and high-mass galaxies (Haines et al. 2007) should be re- 
fiect in variations with mass of the structural and kinemat- 
ical properties and hence in variations both in the shape 
and orientation of the FP for these two families of galax- 
ies. It should also be noted that non-merger origins may be 
important for lower mass galaxies, whose evolution turns 
out to be primarly driven by the mass of their host halo, 
probably through the combined effects of tidal forces and 
ram-pressure stripping (Haines et al. 2007). 

To date no large homogeneous sample of galaxies cover- 
ing both the giant and dwarf regime exists. Although many 
FP datasets for systems as different as BCGs, normal Es, 
dEs, dSphs have been analysed a nd compared to look for 
changes in a (|Zaritskv et al.|[2006l ) , these studies suffer from 
the non-homogenity of the samples both in terms of dif- 
ferences in measuring the quantities entering the FP (for 
example different fits to derive structural parameters, dif- 
ferent apertures to measure the velocity dispersion) as well 
as the selection of galaxies themselves. The form of the FP 
obtained can be infiuenced by all of these criteria. 

In this paper we present a FP analysis of 141 early-type 
galaxies from the Shapley supercluster (2:~0.049) with both 
new i?-band surface photometry measuremen ts and pub- 
lished velocity dispersion measurements from ISmith et al.l 
(|2007l ). The sample is randomly selected down to -I- 3 
and represents the largest homogeneous sample of low-mass 
early-type galaxies with reliable velocity dispersions down to 
~50kms^^. In Section 2, we present the photometric and 
spectroscopic data (including velocity dispersion measure- 
ments and spectral indices) and the catalogue of "newly 
derived" structural parameters; the morphological classifi- 
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Right Ascension (J2000) 

Figure 1. Density map of the SOS area (limited by black lines) with supercluster galaxies observed spectroscopically but not included 
into the FP sample (blue points) and those that enter the FP (red points) indicated. G reen b oxes show the regions covered by images 
of NFPS. The black contours are isodensity contours of R<21 galaxies (sec Haines et al. |2006|) . 



cation is described in Section 3. Our FP fits for the overall 
and high-cr (cr>100kms^^) samples are presented in Section 
4, which discusses how the selection criteria can affect the 
values of the FP coefficients and mimic a possible curvature 
of the plane. In Section 5 we quantify the contribution of 
stellar population to the intrinsic scatter. We discuss our 
results in Section 6 and give a summary in Section 7. The 
origin of the tilt of the FP will be investigated in a forth- 
coming paper. 

Throughout this paper we use Ho=70kms~^ Mpc~^, 
f2m=0.3 and JIa^O.T. With this cosmology 1 
arcsec=0.96 kpc at 2=0.049 and the distance modulus 
is 36.69. 



2 THE DATA 

The sample of galaxies used in this work consists of 141 
early-type _R<18 galaxies (red points in Figure[T} distributed 
thro ughou t entire Shapley Optical Survey (SOS, Mercurio 
et al. |2006| ) area, but mainly located in high-density regions 
(i.e. cluster cores). The sample selection is summarized in 
Table 1. In the region covered by SOS there are 378 con- 
firmed supercluster members. Here we analyse the sample 
of 141 galaxies which have all the following characteristics: 
(i) an early-type morphology (see Section 3), (ii) reliable 
surface photometry, (iii) a measured velocity dispersion and 
(iv) an insignificant Ha emission (EW(Hq)<3A). 

2.1 Spectroscopic data 

Spectra were obtained using the AAOmega fibre-fed spectro- 
graph at th e Anglo-Australian Teles c ope. A full description 
is given by ISmith. Lucev fc HudsonI (|2007l . hereafter SLH) ; 
we summarize some key points here. The spectroscopic sam- 
ple limit is i?=18 in the cluster cores, with galaxies se- 



lected from the N OAO Fundamental Plane Survey (NFPS; 
ISmith et al.l |2004| ) images (green boxes in Figure [TJ . Out- 
side of these regions, brighter galaxies were selected from 
the 2MASS Extended Source Catalogue (i?<15.7). The fi- 
bre diameter of 2arcsec corresponds to 1.9 kpc. Long inte- 
grations resulted in high signal-to-noise ratios, S/N~45 per 
A for cr<100kms"^ and S/N~90 per A for cr>100kms~\ 
Velocity dispersions were measured with respect to the best- 
matching simple stellar population (SSP) templates. The er- 
rors were estimated from Monte Carlo simulations, and are 
~0.05 dex at cr<100 kms"^ and ~0.01 dex for cr>100 kms"^ 
The spectral resolution of 3.2 A (82 km s~^ instrumental dis- 
persion) allows the recovery of velocity dispersions as low as 
~40kms~^. However, some galaxies in the sample remain 
kinematically unresolved, i.e. have velocity dispersions con- 
sistent with zero. When comparing the new velocity disper- 
sion measurements with those previously obtained from the 
NFPS (having a factor ~3 lower S/N), SLH find the new 
velocity dispersions of low-cr objects to be systematically 
lower by ~0.10dex, which they attribute to a combination 
of higher signal-to-noise and the use of a range of SSP tem- 
plates rather than individual K-giant stars. The absorption 
index data are tabulated by SLH. Here we also make use of 
single-burst equivalent ages, metallicities (Z/H) and abun- 
dance ratios (a/Fe) estimated using the models of Thomas, 
Maraston & Bender (2003). Details of this process, and the 
stellar population parameters for each galaxy are provided 
by Smith, Lucey & Hudson (2009). The resulting scaling re- 
lations between the three stellar population parameters and 
velocity dispersion are: 

0.43±0.05 ryirr 0.32±0.04 0.20±0.03 /r,\ 

age oc (T ,Z/iioccr ,a/heoia . (2) 

Typical errors for galaxies with a in the range 50-100 kms~^ 
are 14% in age, 0.05 dex in [Z/H] and 0.04 dex in a/Fe, 
while they reduce to half these values for galaxies with 
(j>150 kms~\ 
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In this work we refer to velocity dispersions measured 
into an aperture of rc/8. We have corrected our velocity 
dispersions ((Tap) acquired with fi bres of 1 arcsec radius (r ap) 
to the apertures of re/8 following I j0re:ensen et all |l993): 

log^ = -0.04 log (3) 

Hereafter, we adopt the notation o"o = o-r^/s for the velocity 
dispersion corrected to this fiducial aperture. 

Starting from the sample of galaxies observed by SLH 
we select those 396 galaxies belonging spectroscopically 
to the Shapley supercluster (0.039<2<0.056, blue and red 
points in Figure 



the images being lower in the overlapping gap area. More- 
over, in the gap area the PSF turns out to be poorly de- 
fined due to its spatial variations. We removed the galaxies 
in these regions from our sample and performed the sur- 
face photometry only for galaxies in highest S/N regions 
(224 galaxies), where the PSF is well defined. We correct 
our mean surface brightnesses for cosmological dimming 
{d{fi)c= 0.208 mag/ arcsec^), galactic extinction [dm— 0.147 
mag, Schlegel et al. 1998) and k-correction {dm = 0.05, 
Poggianti 1997) and convert from magarcsec^^ to log {/)^ 
expressed in physical units L_r0 pc^^ through log{I)^= — 
O.4((^)^-Mfl0-5Iog(2O6265 pc/lOpc)) where Mh0=4.42 is 
the solar absolute magnitude (Binney & Merrifield 1998). 



2.2 Imaging data 

For photometry we refer to the SOS survey (|Mercurio et al.l 
l2006l : lHaines et aLlbOOel ). The SOS is based on data acquired 
with the WFI camera (4x2 mosaic of 2k x 4k CCDs with 
a pixel scale of 0.238" pixels) mounted on the ESO/MPI 
2.2-m telescope at la Silla observatory. 7?-band imaging was 
acquired in good seeing conditions (FWHM~0. 7 arcsec) for 
eight contiguous fields covering a region of 2deg'^ centred 
on the Shapley supercluster core (see Figure [T|. Five expo- 
sures were obtained for each field giving a total exposure 
time of 1200 s for the mosaic images (240 s x 5). The reduc- 
tion was carried out with the ALAMBIC pipeline (version 
1.0, Vandame 2004) an d the catalogue was pro duced with 
the SExtractor package (|Bertin fc Arnoutg|l996l ). plus a set 
of procedures designed ad hoc to remove spurious detections 
(bad pixels, cosmic rays, etc) and correct the photometry for 
blended sources. Th e survey is com plete to R=22 (~M*-|-7). 
For more details see ercurio et al. (2006. ). SOS photometry 
is available for 378 supercluster galaxies observed by SLH, 
all of which are detected at signal-to-noise levels greater than 
100 in each exposure, such that reliable structural parame- 
ters can be derived. 



2.3 Structural parameters 

Structural parameters were derived usi n g the software 
2DPHOT described bv iLa Barbera et all (|2008br ). This is 
an automated tool measuring both integrated and surface 
photometry of galaxies and is furnished with several tasks 
to carry out reliable star-galaxy separation, PSF modelling 
and estimation of catalogue completeness. The main steps of 
the 2DPHOT algorithm are: i) creation of a clean catalogue 
of the input image with SExtractor; ii) estimation of the 
FWHM and the definition of "sure stars"; iii) construction 
of an accurate PSF model taking into account both possi- 
ble spatial variations as well as deviation of stellar isophotes 
from circularity; iv) derivation of structural parameters (ef- 
fective radius r^, mean surface brightness {n}^, Sersic index 
n, total magnitude mtot, etc.) by fitting galaxy images with 
2D PSF-convolved Sersic models, as well as the measure of 
the fit accuracy (xt)- 

The measurement of structural parameters is strictly 
dependent both on the PSF model and on the signal-to- 
noise ratio. Since the mosaic SOS images are obtained by 
stacking jittered images to cover the gap regions between the 
eight CCDs, the signal-to-noise ratio is not constant among 



2.4 Measuring the uncertainties in effective 
radius and mean surface brightness 

To derive the errors, we measure the structural parameters 
on the five single exposures {rl, with i=l, 5), assuming 
that the observed spread of their values is directly related 
to the S/N on the final mosaic used to derive the structural 
parameters entering the FP. 

Due to the size of the sample (five exp osures) only 
two s cale estimators can be considered reliable (|Beers et al.l 
Il990h : the classical standard deviation and the that obtained 
by the gapper algorithm (Jgap. The gapper is a robust scale 
indicator (Wainer & Thissen, 1976) based on the gaps be- 
tween ordered statistics which has a high level of efficiency 
for samples as small as five objects. If we have n measures 
of a quantity x ranked in increasing order (xi,X2,...x„_i,x„), 
according to the gapper algorithm we can measure a robust 
scale indicator as; 



n{n — 1) 



where 



gi = Xi+i — Xi, I ■ 



1 and Wi = i{n — i). 



(4) 



(5) 



To avoid an overestimation of the errors due to the presence 
of outliers, we first compute the o-gap of the log rl and log (7)^ 
distributions and reject all values that deviate more than 
3crgap from the median value, before subsequently computing 
the classical standard deviation as well as the covariance 
matrices of the clipped sample for both variables (cj^^g , 
^iog(/>j. cov[logro;log{7)J). 

The errors on the mosaic value of logrc and log (7)^ 
(51ogrc and 5iog(/>„ respectively) are given by a^^^^i j yfn 



and CTiQg i^jy I 'Jn^ where n is the number of measures avail- 
able. The typical errors on the logarithms of effective radius 
(logrc) and mean surface brightness (log(7)j,) are 0.03 and 
0.04, respectively. These errors explicitly include the effects 
of noise in the galaxy, but not the presence of neighbour- 
ing objects. We have estimated the effect of the latter, by 
repeatedly placing copies of the galaxies one-by-one at ran- 
dom positions across the same CCD image (where the PSF 
should remain constant) and reapplying 2DPH0T, finding 
the variations in the structural parameters to be consistent 
with the previously obtained errors, albeit with a small num- 
ber (~2%) of >5(T outliers when the galaxy is placed very 
close to a bright star or galaxy. 
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Figure 2. Upper panels: Images of six early-type galaxies covering the range of magnitudes of our sample with their corresponding 
residual maps beneath. Lower panels: equivalent for six late-type galaxies. The residual maps facilitate detection of spiral arms and other 
irregularities in the late-type ojects. These examples demonstrate the success of 2DPHOT in fitting the galaxy profile even when the 
object is near the edge of the image and/or bright objects. 



3 MORPHOLOGICAL CLASSIFICATION 

We have morphologically classified the 224 galaxies with 
available surface photometry, by inspection of the residual 
maps provided by 2DPHOT. We denote as "late-type" all 
galaxies showing signs of spiral arms or asymmetric distur- 
bance, and as "early-type" all those with no such structures. 
The resolution of SOS images does not allow any finer clas- 
sification, for example into elliptical and lenticular galax- 
ies, since the presence of a residual disk can be seen only 
in particular cases, i.e. when it is very bright, widespread 
or edge-on. In Figure [2] we show some illustrative exam- 
ples of galaxies classified as early- and late- types, covering 
the full range of magnitudes studied here. We have checked 
this classification by com paring our results with those of 
iThomas fc KatgertI (|2006t ) for a subsample of 54 galaxies 
belonging to A3558 and A3562, finding perfect agreement. 
In their classification galaxies were subdivided into E, SO, 
Se and SI classes (the latter being early and late spirals). 
Our early- type sample contain only galaxies classified as E 
and SO bv lThomas fc KatgertI (|2006l ). 

The str ong correlation betwee n the Sersic index n and 
luminosity (|Young fc Currid Il994l : ICaon et al.1 1 19931 ') pre- 
vents any morphological classificatio n based on Sersic in- 
dex alone. iGraham fc GuzmanI (|2003l ) show that the values 
of Sersic index n of dE, ordinary E/SO galaxies and BCGs 
follow a continuous trend from 7i<0.5 to 7i~10, whereby 




Figure 3. Distribution of Sersic index as function of morpholog- 
ical type and magnitude. The filled circles indicate the median 
value of n for early-type galaxies in each magnitude bin (chosen 
to contain the same number of galaxies), while the "error bars" 
indicate the interquartile values. The open circles correspond in- 
stead to late-type galaxies. 



brighter galaxies have larger values of n. In Figure |3] we 
demonstrate the equivalent trend for our galaxy sample, 
plotting the Sersic index distribution as function of mor- 
phological type and magnitude. The median value of n for 
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Table 1. Spectroscopic sample: the first four fines indicate the cri- 
teria used to select the galaxy sample, and the number of galaxies 
remaining after each step. In the foUwoing lines, the 224 remain- 
ing galaxies are then divided into four classes, of which only the 
141 early-type objects are used in our FP analysis. 



Sample selection 


Galaxies with spectra 


565 


Cluster members {0.039<z<0.055) 


396 


Cluster members in SOS region 


378 


Cluster members in SOS region with surface photometry 


224 


Galaxy classification 


early-type 


141 


late- type 


44 


unresolved kinematically 


31 


with EW(Ha)>3A 


8 



early-type (filled circles) and late-type (open circles) galax- 
ies are shown for each magnitude bin (chosen to contain the 
same number of galaxies) , while the "er ror bars" indicate the 
inter- quart ile values. As pointed out bv lGraham fc GuzmanI 
the distribution of n for early-type galaxies is shifted 
towards higher values with increasing luminosity. Late-type 
galaxies have lower values of n at all magnitudes, but there 
is not a clean a priori separation between the two morpho- 
logical classes. 

The morphological classification results in a sample of 
141 early-type galaxies and 44 late-type or irregular galax- 
ies, while 31 galaxies are removed from the sample since they 
are not kinematically resolved, and 8 are removed as having 
significant Ha emission (EW(Hq)>3A). We summarize the 
characteristics of the sample in Table [T] while in Table [2] we 
report structural and kinematical parameters with the cor- 
responding errors for the 141 galaxies used in the following 
section. 



4 THE FUNDAMENTAL PLANE OF R<18 
SHAPLEY GALAXIES 

For the FP we use the representation of Equation [T] where 
Te is the effective radius measured in kpc, ao is the central 
velocity dispersion corrected to an aperture of radius ro/8 
in kms^^ (see Section 2.1) and (I)^ is the mean surface 
brightness within To expressed in Lq pc~^ . To derive the 
value of a, /3 and 7 we adopted the orthogonal fit which 
minimizes the quantity: 



E 



|logr^-(alogag-f /31og (J): + 7)! 
^1 + q2 -h /32 



(6) 



i.e. the sum of the absolute residuals perpendicular to the 
plane. This is les s sensitive to outliers than the classic least- 
squares method (jjeirgensen. Franx fc Kiaegaardlll996l . here- 
after JFK96). For the full sample of 141 galaxies we found: 



logr-e=1.03±0.06 log (7o-0.82±0.02 log (1)^+0.33, 



(7) 



where the errors on coefficient are computed via a bootstrap 
procedure. The above relation presents a scatter of 0.088 dex 
in the log r^ direction. We estimate the intrinsic contribution 
to the scatter by subtracting in quadrature from the rms 
scatter, the observational errors on each variable taking into 
account the correlation between the errors on the effective 
radius and surface brightness. We obtain an intrinsic scatter 
fint = 0.068. When deriving the (7i„t the errors on distance 
estimated through the standard deviation of the redshift 
distribution are also considered. 

To take into account the different uncertainties on ctq, 
re, along the plane , a weighted fit is adopted. Following 
ICappellari et al.l (|2006l ) we add in quadrature to the mea- 
surement errors the intrinsic scatter orthogonal to the FP, 
5int±. The fitting function to be minimized is: 



^ Wi\ log rl-a log ao- 13 log {1)1 - 7! , 



(8) 



where the sum is extended over all sample galaxies and 
where 



-'-Si 



-(^logr, +(a5iog<j(,) +(/35iog/J -2/3cov[logro;log(J)J.(9) 



The intrinsic scatter, 5int±, is adjusted to give the expected 
value of 0.8 per degree of freedom (the mean absolute value 
of a standardised gaussian distribution is 0.8). The resulting 
equation for the FP is: 



log r-o=1.06±0.06 log o-o-0.82±0.02 log (1)^+0.28. 



(10) 



The orthogonal and weighted fits are consistent, demon- 
strating that the results are robust to the effects of the larger 
ao uncertainties (0.01 dex for ao > 100 km s~^ and 0.05 for 
ao < 100 km s~^) for the lowest mass galaxies. The scatter 
in the log re direction is now equal to 0.092, while its in- 
trinsic component is 0.070. Figure |4] shows the two edge-on 
FP views (left and central panels) and its trend along the 
direction of luminosity (right panel) . 

The value of the a coefficient obtained for the Shapley 
sample, is significantly lower than the typical values of 1.2- 
1.3 reported for samples dominated by giant galaxies (e.g. 
JFK96). If we restrict our sample to galaxies with ao > 
100 km s"'^ (high-(Jo sample, 91 galaxies) we obtain the FP 
relation: 



logr-o=1.35±0.111og ao-0.81±0.03 log {1)^-0.40, 



(11) 



with the overall scatter in logr-o reduced to 0.090. In Fig- 
ure [5] we plot the FP as fitted for the high-ao sample (black 
points). Galaxies with ao<WO kms^^ are also reported for 
comparison (red dots). The a value of the high-cro sample is 
closer to those of other authors. The lower a value found for 
the total sample is thus probably due to the extension of our 
sample to very low mass galaxies, down to ao ~50kms~^. 
The improved method used to obtain the velocity disper- 
sions, as well as the high signal-to-noise levels of the spec- 
tra, can produce systematic effects on the resultant values 
of ao, particularly for low-ao objects. SLH find that velocity 
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Table 2. Structural and kinematical parameters for the FP galaxy sample. Cofamn 1: ID (SLH); Columns 2 and 3: RA(J2000) and 
DEC(J2000); Column 4-'- total i?,-band magnitude; Columns 5 and 6: log ctq referred to an aperture of re/8 radius and i5iogo-o; Columns 



7 and 8: log rc and i5|ogro) Columns 9 and 10: log (I)^ 
be available in the electronic version. 



and <5i 



log (/> ; Column 11: n; Column 12: reduced xt- The whole catalogue will 



ID 


RA(.T2000) 


Dcc(.72000) 


R-tot 


log o-f, 




log 








n 












(km/s) 




(kpc) 














TMASSJ13274662-3059237 


13:27:46.6 


-30:59:24 


15.166 


2.201 


0.010 


0.904 


0.042 


1.766 


0.070 


10.708 


1.25 


NFPJ132828. 6-313205 


13:28:28.6 


-31:32:04 


15.404 


2.112 


0.012 


0.308 


0.002 


2.861 


0.004 


1.589 


6.43 


NFPJ132738. 0-313041 


13:27:38.0 


-31:30:41 


16.549 


2.044 


0.016 


0.038 


0.004 


2.943 


0.008 


2.283 


1.31 


NFP.I133408. 2-314735 


13:34:08.1 


-31:47:34 


16.083 


1.757 


0.097 


0.682 


0.044 


1.842 


0.059 


3.518 


2.03 


TMASSJ13295423-3201001 


13:29:54.2 


-32:00:59 


15.229 


2.153 


0.010 


0.475 


0.005 


2.597 


0.007 


4.617 


1.17 


NFPJ132656. 0-312528 


13:26:56.0 


-31:25:27 


14.813 


2.031 


0.007 


0.797 


0.007 


2.118 


0.011 


5.452 


1.99 


TMASSJ13294834-3115580 


13:29:48.3 


-31:15:58 


14.075 


2.296 


0.004 


0.984 


0.023 


2.043 


0.035 


6.945 


1.15 


NFPJ132810. 5-312310 


13:28:10.5 


-31:23:09 


13.868 


2.267 


0.004 


1.108 


0.043 


1.875 


0.067 


7.796 


1.48 


NFPJ132426. 5-315153 


13:24:26.5 


-31:51:53 


14.857 


2.302 


0.006 


0.506 


0.036 


2.685 


0.051 


8.414 


1.38 


NFPJ132923. 8-314832 


13:29:23.2 


-31:48:32 


15.060 


1.718 


0.062 


0.972 


0.046 


1.673 


0.073 


9.545 


1.77 



C7> in 

5 d 



' ' ' I ' ' ' ' I ' ' ' ' I ' 

a= 1.06 +/- 0.06 

|3 0.82 +/- 0.02 

y= 0.28 +/- 0.12 




Total sample (R<18.0) 
-r 



0.5 1 1.5 2 
1.06 Log (7o - 0.82 Log <\>^ + 0.28 





18 16 14 12 

(a/^) Log <7o - (1+2/!)/|8 Log r. + (r//S) + k 



Figure 4. Left, and central panels: The edge-on views of the Shapley i?<18 FP. Right panel: The edge-on view of the FP as it appears 
along the direction of luminosity. Black lines are the best-fit relations, a, (3 and 7 values are reported in the left panel. 



dispersions obtained using single, old, solar-metallicity mod- 
els instead of templates with different metallicities (as used 
in this paper) are overestimated by ~6% for (Jo=75kms~^ 
galajcies and by ~18% for (To=50kms^^. As a result, in our 
sample the low-cro values are systematically lower than those 
of previous samples (e.g. NFPS sample) , causing an increase 
of the FP tilt of about 10-15%. The impact of the Iow-cto 
limits is explored further in Section 4.2. We notice that re- 
sults presented in Section 5 and discussed in Section 6 has 
been checked to be independent of ctq cuts. 



4.1 Comparison with Coma 

Besides the range of velocity dispersions, the determina- 
tion of the three FP coefficients is strictly dependent on 
other factors such as the selection criteria of the sample, the 
fit algorithm and the proced ure used to derive structural 
and kinematical parameters (|Kelson et al.l I2OO0I ) . We ad- 
dress this problem comparing our FP with the previous work 
by JFK96. Their sample consists of 81 early-type galaxies 
down to Gunn r<15.3 in the central region of Coma clus- 
ter. For this section, we have corrected our velocity disper- 
sions to their standard fixed aperture of 0.595 kpc of 
radius. To analyse the consistency of the two samples, in 



Figure |S] we plot the spectroscopic completeness functions 
(CFs) of the Coma (open circles) and Shapley (filled cir- 
cles) FP samples as a function of magnitude. The magni- 
tudes of the Shapley galaxies are converted from R to the 
Gunn-r photometric system accordin g to the typical colou r 
of r— i?=0.35 for elliptical galaxies (|Fukugita et al.lll995l ). 
The Coma CF is computed using the JFK96 FP sample 
and the complete pho tometric catalogue of early -type Coma 
galaxies published bv lJorgensen fc Franxl (|l994h . The limit- 
ing magnitude of the Coma sample is Mr=— 19.79 (dashed 
line). The Coma sample is complete at bright magnitudes, 
but the completeness declines rapidly towards zero at the 
faint magnitude limit. In contrast, about 30% of Shapley 
galaxies are spectroscopically observed and have available 
velocity dispersion measurements independent of magnitude 
down to Mr=— 19.79. In the magnitude bins fainter than 
Mr = —20, the Shapley CF declines gently due to an in- 
creasing fraction of galaxies with surface brightnesses and 
velocity dispersions too low for successful (To measurements, 
before dropping rapidly in the faintest bin. Green triangles 
represent the CF of Shapley galaxies observed spectroscop- 
ically with no regards to successful velocity dispersion mea- 
surements: the decline at faintest magnitude has now disap- 
peared except for the faintest bin (Mr>— 18.5). The distri- 
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Shapley: > 100 km/s sample (R<18.0) 




1.35 Log <7o - 0.81 Log <\>^ - 0.39 (Log - |8 Log <\>^ - y)/a {a/fi) Log - (1+2jS)/|8 Log + (yZ/S) + k 

Figure 5. Left, and central panels: The edge-on views of tlie Siiapley R<18, (to>100 km (black dots) FP. Right panel: The edge-on 
view of tlie FP as it appears along tlie direction of luminosity. Black lines are the best-fit line of the high cr sample. Galaxies with 
o-o<100kms^^ (red dots) are reported for comparison, a, 13 and 7 values are reported in the left panel. 



T ' r 




Mr 

Figure 6. The spectroscopic completeness function both for 
Coma (open circles) and for Shapley (filled circles) sample. Green 
triangles show the completeness for all Shapley galaxies observed 
spectroscopically, including those with no successful measurement 
of cr. 



Shapley Mr < —19.79 Log cJq > 2.0 galaxies 




Figure 7. Edge on view of the FP for Shapley early-type galaxies 
with Mr.< — 19.79 and (To>100kms~^ (Shapley "matched" sam- 
ple, black points). The solid line represents the Shapley FP fit. For 
comparison are reported also Coma galaxies from JFK96 (open 
points). 



bution reflects well the random selection criteria of our spec- 
troscopic survey for R<18 galaxies. To compare the FPs of 
the two samples we have selected from Shapley only early- 
type galaxies with Mr< — 19.79 and (7o>100kms~^ (here- 
after "matched" sample) corresponding roughly to the lim- 
its of the Coma FP sample. This new sample consists of 88 
galaxies. The FP of the "matched" sample is: 

logr-o = 1.21±0.081ogcro - 0.75±0.02 log (J)^ - 0.19, (12) 

with a scatter equal to 0.08 dex in log rc direction. In Fig- 
ure [7] we plot the edge-on view of the FP for the "matched" 
sample (filled circles). Coma galaxies (open circles) are re- 
ported for comparison. Since in JFK96 the value of Ho is 
set equal to 50kms~^ Mpc~^, the effective radii of Coma 
galaxies have been shifted by a factor -0.146 dex. Both the 
distributions and the dispersions of the two samples are con- 



sistent, in fact the FP for Coma sample, as found by JFK96 
is: 

log r-o=1.31±0.07 log (7o-0.84±0.02 log (1)^-0.082. (13) 

The possible reasons for the remaining slight differences be- 
tween the two FPs could be the different sampling at the 
brightest magnitudes of the two datasets, and the different 
procedures used to obtain the structural parameters, since 
JFK96 fit their early-type galaxies with de Vaucouleurs pro- 
files. 

4.2 A curved surface or a selection effect? 

Observing the high-ao fit (see Figure [5]) it is notable that 
all galaxies with velocity dispersions less than lOOkms"^ 
(low-mass galaxies) are systematically displaced above (left 
panel) or below (central panel) the best-fit plane as defined 
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CTq cut off (km/s) 



Figure 8. Values of the a coefficient for Sliapley sample of galax- 
ies with increasing lower cuts on velocity dispersion. The point 
with (To<50kms~^ represents the coefficient for the whole sam- 
ple. 



by the high-cro sample. In particular, the log cto edge-on view 
suggests that low-mass galaxies do not follow the same rela- 
tion as massive galaxies. Thus, at face value, the FP appears 
curve d, as s uggested also by other studies (JFK96, Zaritsky 
et al. I2OO6I . D'Onofrio 2008, Nigoche-Nietro 2008) investi- 
gating the curvature of the FP as a function of mass and/or 
luminosity. Studying a sample of 69 faint early-type galax - 
ies in the core of Coma cluster. [Matkovic fc GuzmanI (|2005l ) 
find a similar behaviour for the L — a relation, with dwarf 
galaxies having a trend (1/ oc a^ '^^**^ '"') shallower than that 
observed for elliptical giant systems (L oc cr^) . 

In Figure [S] we investigate this possible curvature by 
analysing the dependence of a as a function of different low- 
(To limits. We cut the galaxy sample at various values of 
(To and estimate the FP coefficients through the weighted 
orthogonal fit procedure. One can see that a increases sys- 
tematically as the (To cut moves to higher values, remov- 
ing more galaxies from the sample. The values of a for the 
whole sample (q:=1.06±0.06) and for the highest-iro cuts 
(a=1.72±0.28) differ by more than twice the standard er- 
rors, indicating a possible curvature of FP. To investigate if 
these observed variations of a with the sample selection are 
due to a real curvature of the FP, or induced by selection 
effects, a set of simulations were performed. We constructed 
1000 mock catalogues using the observed values of log (7)^, 
and log (To and the value of log assigned using the relation: 



log r-o = a log (TO + /3 log (J)^ + 7 + ^^(0, S), 



(14) 



where, in this case, a = 0.92, /3 = -0.78, and 7 = 0.496 
are the coefficients obtained by fitting the FP of the over- 
all Shapley sample with log as the dependent variable, 5 
= 0.099 is the observed scatter in the log re direction, and 
N(0,(5) is a Gaussian random variable with zero mean and 
standard deviation of 5. For each mock catalogue we evalu- 
ate the relative change in a {5a/ a) between the orthogonal 
fits obtained for the whole sample and after applying a cut 
at (To = lOOkms"^. The resulting histogram of 5a/ a. val- 
ues is shown in Figure [51 where the dashed line indicates 
the 5a/ a computed directly for the Shapley sample. This 
shows that the observed change in a is fully consistent with 




5a/a 

Figure 9. The distribution of the relative variation of a coeffi- 
cient after applying the cro > lOOkms"^ cut for the simulated 
sample. The dashed line marks the variation of a when applying 
the (TO cut to the Shapley data. 



o 
b 



LD 
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o 




Log - /? Log <l>^ 



7 



Figure 10. The edge-on projection of the FP for the real Shapley 
sample (black dots) and for one of the simulated catalogues (open 
circles). The line represents the best-fitting relation for the subset 
with cro > 100 km s"'^. 



that expected for a linear relation. Moreover, in Figure [TO] 
we compare the FP edge-on projection (Eq. 11) for the ob- 
served high-(To sample (solid symbols) and a simulated sam- 
ple obeying the overall Shapley FP relation (open symbols). 
The two samples show similar behaviour, with the same ap- 
parent curvature such that galaxies with low velocity disper- 
sions are placed systematically below the FP relation. Given 
that the simulated sample is defined to follow a linear rela- 



10 Gargiulo et al. 



tion, this suggests that the apparent curvature in Figure [5] 
can be explained solely by selection biases. Hence, in spite of 
the wide-range in velocity dispersions covered by our sample, 
this is not sufRcient to distinguish between either a linear or 
curved FP relation, but there is no convincing evidence for 
the latter in our data. 



5 FP RESIDUALS AND STELLAR 
POPULATIONS 

We find that for the Shapley FP, the contribution of the 
intrinsic scatter (0.070) to the total scatter (0.092) is larger 
than that of the measurement errors (0.060). Taking ad- 
vantage of the wide variety of indices derived by SLH we 
investigate the possible origin of the intrinsic scatter of the 
FP analysing the correlations between both spectral indices 
(Mgb, Fe5015, Fe4383, H/3, HgF, HdF) and stellar popula- 
tion parameters (age, metallicity, a-enhancement) and the 
residuals from the FP. 



5.1 FP residuals vs. single spectral indices 

All of the many spectral indices are known to correlate with 
(To, which, if not corrected for, could produce spurious corre- 
lations with the FP residuals. To avoid this problem, rather 
than use the spectral indices themselves, we consider their 
residuals with respect to the index-log cro relation. For each 
index, we fitted first the index-log ao relation, shown in Fig- 
ures [11] and [12] (bottom panels) for metallicity- and age- 
sensitive indices respectively, and then determine the resid- 
uals with respect to that relation (middle panels). In the 
figures, the upper panels show the index-log ao residuals 
versus the residuals from the FP in the log re direction. For 
each relation, we quantify the product-moment correlations 
between the residuals and spectral index through the corre- 
lation factor r and the bisector least-squares fit, assuming 
that the distributions are both Gaussian. Uncertainties in r 
and the fits are estimated through 10000 Monte Carlo real- 
izations, taking into account the errors on each value. These 
values together with the probability p that two quantities 
with correlation factor r are not correlated are reported in 
Figures [Hand [H 

Among the primarily metal-sensitive indices (Figure 
11), the strongest dependence is found for Fe4383, which is 
correlated with FP residuals at the 4a level. A similar but 
weaker correlation is obtained for Fe5015. For Mgb, how- 
ever, a weak correlation is obtained in the opposite sense. 
Thus galaxies which are more compact than expected from 
the FP have lower Fe4383 indices than expected for their ao, 
but higher Mgb indices than expected for their ctq. This al- 
ready provides a hint that the physical driver for the residual 
trends is not the total metallicity, but instead the ratio be- 
tween Mg and Fe abundances. The Balmer lines show weak 
positive correlations, such that galaxies with strong index 
values (for their ao) are more diffuse than predicted by the 
FP. This could be caused either by age or by metallicity ef- 
fects, since the Balmer lines are sensitive to both parameters 
to some extent. 



5.2 FP residuals vs. stellar population parameters 

Using a single spectral index it is not possible to disentangle 
between the effects of age and metallicity. In this section, we 
use estimates of age, metallicity and a-enhancement (a/Fe) 
derived from the index measurements by means of a multi- 
index procedure, to provide more physically meaningful in- 
formation (Smith et al. 2009). 

As for the single spectral indices, the three stellar 
population parameters age, metallicity and a/Fe correlate 
strongly with ao as shown by the scaling relations of Equa- 
tion [21 therefore we first compute the residuals of stellar 
population parameters with respect to these relations. The 
lower panels of Figure [13] show the stellar population residu- 
als against log cro characterized, by construction, by the ab- 
sence of correlations between the residuals and log ao itself. 
The top panels show the resultant correlations between the 
age, metallicity and a/Fe residuals (left, middle and right 
panel respectively) and the FP residuals along the log rc 
projections. 

Both age and a/Fe are seen to be strongly anti- 
correlated (at >3a and >4(7 respectively) with the residuals 
computed in the log rc direction, while metallicity shows a 
positive correlation with the FP residuals along the same 
axis. We find these correlations to be independent of the 
range of a covered, for example obtaining consistent results 
when using only the high-ao sample. In Table [3] we summa- 
rize the probability p for all the cases. 

We note that in the figures 11, 12 and 13, the veloc- 
ity dispersion enters into both axes, since each quantity is a 
residual from a relation involving ao- The errors are hence 
correlated which will produce a bias towards positive cor- 
relations. For each of the primary SSP parameters we mea- 
sure the level of this bias through Monte Carlo simulations, 
based on fake data in which there are no intrinsic corre- 
lations. Firstly, for each galaxy we assign fake SSP values 
(age, Z/H, a/Fe) following the previously derived index-ao 
relations, and having the same rms intrinsic scatter in the 
SSP value about the relation. We then perturb each of the 
values log cro, log rc, and log (/c) by their uncertainties for 
each galaxy, and recalculate the residuals about the FP and 
the index-cro relation as before. Through these simulations, 
we measure this bias (A r), i.e. the difference in the corre- 
lation coefficient, due to the correlated measurement errors 
in cro and finds: 

Ar(logage) = 0.03, Ar(Z/H) = 0.02, Ar(a/Fe) = 0.04.(15) 

We see that this bias contributes to the observed correlation 
between the metallicty and FP residuals, but acts to op- 
pose the anti-correlations seen for age and a-enhancement. 
In each case the effect of the bias is small, being at a level 
of ~10% of the observed correlation. The relation between 
residuals obtained taking into account this bias is shown in 
figure 13 by the red dot-dashed line. 

It should be remembered that the stellar population pa- 
rameters presented here are derived through fixed apertures 
covering only the galaxy centres. Therefore, if galaxies have 
significant population gradients, the observed trend could 
be due to sampling the central stellar populations at differ- 
ent radii. In the case of metallicity, this may be a factor, 
as early-type galaxies are observed to have negative metal- 
licity gradients (Kuntschner et al. 2006; Rawle et al. 2008) 
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Figure 11. Low-left panel: The Fe 4383-logo'o relation. In all the panels, the solid (dashed) lines indicate the mean (and Icr confidence 
limits) bisector least-squared fits after 10000 Monte Carlo realizations. The corresponding correlation coefiBcient r (and uncertainty) 
is indicated. Middle-left panel: The corresponding Fe 4383-logo'o residuals in Fe 4383 direction respect to the fit (solid line in low-left 
panel) plotted against logcro. Top-left panel: Correlations between the residuals of the Fe 4383- log (tq relation in Fe 4383 direction with 
the residuals from the best-fitting FP in the log r^ direction. The probability p that two quantities with correlation factor r are not 
correlated is also reported. Central panels: The same for Fe 5015. Right panels: The same for Mgb. 



Table 3. Probability that two quantities with correlation factor r as in Figure 11, 12 and 13, are not correlated. 



Mgb 


Fg5015 


Fc4383 


Hbcta 


HgF 


HdF 


Log ago 


Z/H 


a/Fc 


11.3% 


5.3% 


0.092% 


2.3% 


14.9% 


4.0% 


0.056% 


0.38% 


0.001% 



of the order —0.20 ± 0.05 dex in [Z/H]. However, spatially- 
resolved spectroscopy of early-type galaxies show that they 
generally have flat radial trends in age and a/Fe (Sanchez- 
Blazquez et al. 2007; Rawle et al. 2008) , indicating that the 
anti-correlation of age and abundance ratios with the FP 
residuals in the log rc and log (7)^ directions is robust. 

The strong correlations found between the stellar pop- 
ulation parameters, principally the a/Fe, and the residuals 
of the FP suggest that its scatter is in part due to variations 
in the stellar populations at fixed galaxy ao- To this aim we 
fit the modified FP relation: 



log r-e = a log ao + P log {I)^ + Kspp^ + 7, (16) 

where the stellar population parameter sppi is, in turn, a 
enhancement, age and metallicity of the galaxies, and study 
the variation induced in the scatter by keeping a and (3 
fixed as in Eq. [14] and allowing k to vary to minimize the 
FP scatter. It is necessary to keep a and j3 fixed as each of 
the stellar population parameters correlate strongly with ao, 
and hence would introduce spurious variations in a. The re- 
sults obtained are listed in TablelH where the second column 
indicates the estimated strength of the spp term k, the fol- 
lowing two columns indicate the overall scatter in the log Ve 
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Figure 12. As for Figure 11, but now showing correlations for the Balmcr Une indices HdF (left panels), HgF (central panels) and H/3 
(right panels). 

Table 4. The impact on the overall and intrinsic FP scatter of adding a further term (ssp) to the orthogonal fit of the fundamental 
plane. Column 1: The k value which minimizes the overall FP scatter along the log direction to Eq. 1161 Column 2: total scatters 
around the FP; Column 3: estimated intrinsic FP scatter; Column 4 estimated rms contribution of the spp term in the log re direction. 



spp 




rms in log ro dir. 


intrinsic FP rms 


spp rms contribution 






0.088 


0.068 




a/Fe 


+0.582 ±0.128 


0.075 


0.049 


0.047 


Log age 


+0.169 ±0.057 


0.084 


0.063 


0.026 


Z/H 


-0.160 ±0.061 


0.087 


0.066 


0.016 



direction after the addition of the spp term, and the relative 
intrinsic scatter, and the final column the rms contribution 
of the ssp term to the scatter in log re. The fit with the 
most significant k coefficient (4.5ct) is the fit with the a/Fe 
as fourth parameter. In the case of age and metallicity, the 
new coefficient is non-zero at just the ~3cr level, indicating 



the weaker impact of these parameters on the fundamental 
relation. 

While it seems that none of the additional spp terms 
reduce the overall scatter significantly, the relative impact 
of the three spp terms becomes clearer when considering 
the intrinsic scatter (i.e. after accounting for the measure- 
ment uncertainties) . Both the relation including the effect of 
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Figure 13. As for Figure 11, but now showing correlations for the SSP parameters, age (left panels), total mctallicity (central panels) 
and a/Fe (right panels). The red dot-dashed line shows the relation obtained taking into account the correlated errors (sec text). 




Figure 14. Illustration of the a/Fe dependence of the FP residuals. In the left-hand panel we define three subsamples using residuals 
from the a/Fe-loga relation. The 50% of galaxies closest to the relation are plotted in grey, while the 25% with a/Fc higher (lower) 
than expected for their a are plotted in red (blue). The central panel shows the edge-on projection of the FP colour-coded according 
to this scheme. It is readily seen that high- A [a/Fe] galaxies are offset from the mean FP in a direction which may be interpreted as 
smaller effective radius, lower surface brightness, or higher velocity dispersion. The low- A a/Fc galaxies are offset from the mean FP in 
the opposite direction. In the right-hand panel, we show the FP after correcting for the a/Fe correlation. A modest reduction in total 
scatter is seen, despite the introduction of uncertainty due to a/Fc measurement errors. 
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age and metallicity have an intrinsic scatter of ~0.06, only 
marginally lower than the value without the spp term, while 
in the case of a/Fe the intrinsic scatter is reduced to 0.049. 
This is comparable to the rms contribution from the a/Fe 
term (0.047), i.e. a/Fe contributes around half of the intrin- 
sic scatter, and indicates that the distribution of galaxies 
around the FP are tightly related to the enrichment, and 
hence to the timescale of star-formation. The a/Fe depen- 
dence of the FP scatter is illustrated further in Figure [TJ] 
which demonstrates also the reduction in total scatter de- 
spite the introduction of measurement errors in a/Fe. We 
split the Shapley sample into three groups (shown as red, 
grey and blue points in Figure I14|l according to the posi- 
tion of each galaxy with respect to the a/Fe-logao relation 
(see left panel): the 50% of galaxies closest to the relation 
are plotted in grey, while the 25% with a/Fe higher (lower) 
than expected for their a are plotted in red (blue). A modest 
reduction in total scatter is seen, despite the introduction of 
uncertainty due to a/Fe measurement errors. 

Furthermore the total FP scatter is reduced by includ- 
ing the extra terms, in spite of the fact that the SSP pa- 
rameters are subject to their own substantial uncertainties. 
This implies there must be a still greater reduction in the 
intrinsic FP scatter. 



6 THE ORIGIN OF THE INTRINSIC 
SCATTER 

JFK96 pointed out that the dispersion around the FP rela- 
tion is not completely due to the measurement errors but has 
an intrinsic scatter whose nature is not yet understood. The 
existence of this intrinsic scatter was interpreted as due to 
another "fundamental" parameter characterising the fam- 
ily of early-type galaxies. In many of the previous works 
on the FP (see Section 1), a strong effort has been made 
to find correlations between the FP residuals and different 
line-indices considered to be representative of a particular 
stellar population parameter (for example H/3 for age, Mg2 
for metallicity etc). 

In Section 5 we showed that the FP residuals are 
strongly correlated with estimated stellar population param- 
eters especially with the a-abundance ratio (a/Fe). Given 
that age and a/Fe are known to correlate strongly for galax- 
ies of a given mass, while age and metallicity anti-correlate 
(e.g. Proctor & Sansom 2002; Smith et al. 2007b), it is not 
clear which of the three parameters drives the correlations. 
However, the fact that the strongest correlation is observed 
for the a/Fe suggests that this parameter is playing the ma- 
jor role, while age and metallicity trends just reflect their 
mutual correlations. The trends observed for the metallicity 
indices also suggest that the a/Fe is driving the correlations. 
In fact, the positive correlations between the Z/H dependent 
indices Fe4383 and Fe5015 combined with the negative cor- 
relation for Mgb, is consistent with an anti-correlation for 
the abundance ratio a/Fe. Moreover, we found that about 
half of the FP intrinsic variance is due to variations in 
the stellar populations as described by the a-enhancement 
trend, while the other stellar population parameters seem to 
contribute rather less to the intrinsic scatter. 

It has long been realised (Tinsley 1980) that in the 
study of galaxy formation and evolution, a major role can 



be played by the analysis of the abundance ratios, due to 
them being relatively model independent and primarily af- 
fected by stellar nucleosynthesis and the initial mass func- 
tion (Matteucci 1996). In particular, the ratio between the 
so-called a elements (synthesized in type II supernovae) 
and iron (mainly from delayed type la events) is widely ac- 
cepted to be affected by the characteristic timescale of star- 
formation, the a elements and iron having different produc- 
tion timescales (Greggio & Renzini 1983). According to the 
simplest and widely accepted scenario, a galaxy with a high 
value of a/Fe has experienced many type II SNe events, but 
almost no SNe of type la during the major epoch of star- 
formation: this constrains the timescale for this epoch of 
stellar formation to be shorter than 3x10* years. 

Our results suggest that the galaxies which formed 
their stars over a shorter duration (high a/Fe) are also 
those which are more compact. Such a pattern may be ex- 
pected in a dissipational merger inducing a nuclear stur- 
burst. Vazdekis, Trujillo and Yamada (2004), studing 21 
early-type galaxies already claimed the existence of a corre- 
lation between Sersic index and Mg/Fe ratio (i.e. between 
galaxy structure and stellar population), even if they did 
not correct for correlations with a. They interpreted their 
results as more massive galaxies having their star formation 
quenched on shorter time-scales. Although, in agreement 
with our finding, their analysis is based on to a smaller sam- 
ple spanning a narrower range of velocity disperisons and 
luminosities with compared to our Shapley sample making 
any comparisons d i fficult . 

iHopkins et al.l (|2008l 'l make specific predictions for the 
effects of varying the dissipational fraction on the remnant 
spheroids of a given mass, resulting in correlations between 
the structural and stellar population parameters. By con- 
sidering identical progenitor disks (at t = Q) with initial gas 
fractions fgaa ~ 1 following an exponential star-formation 
history with time-scale r, we have that the gas fraction at 
the time of the merger tm (and hence dissipational frac- 
tion in the merger-induced starburst) will scale as /gas = 
exp(— tm/T"). If the remaining gas is then consumed in the 
central star- burst, producing a passively-evolving spheroid 
remnant, then the dissipational fraction wi ll increase for 
earlier mergers, and hence mean stellar age. iHopk ins et al.l 
(|2008l ) also show that ellipticals formed through mergers 
with higher dissipation fractions should be more a-enriched. 
As merger remnants involving more dissipation, like those 
between more gas-rich disks, are expected to yield larger 
mass fractions formed in nuclear starbursts, which reduce 
significantly the effective radii of the remnant, we should ex- 
pect both mean stellar age and a/Fe to anti-correlate with 
the residuals in the direction. 

We observe clear correlations between the FP residuals 
and age and a/Fe, with those galaxies of a given mass with 
effective radii smaller than predicted by the FP to have stel- 
lar populations systematically older and with higher abun- 
dances than average, fully consistent with the predictions of 
iHopkins et all (|2008l ). 



7 SUMMARY AND CONCLUSIONS 

We have derived the FP of a sample of 141 early-type i?<18 
galaxies in the Shapley supercluster at z=0.049. Velocity 
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dispersions and stellar population parameters were derived 
from the spectroscopic data of Smith et al. (2007), while 
R-band photometry is from the Shapley Optical Survey 
jMercurio et all 120061 ). The final sample extends down to 
in magnitude and 50 km s^^ in ctq. Using the soft- 
ware 2DPHOT (La Barbera et al. 2008b) we derived for each 
gala^xy the structural parameters ro, (fJ-)^, n, and mtot by fit- 
ting a 2D PSF-convolved Sersic model. The morphological 
cl assification was perf o rmed by eye and checked with those 
of lThomas fc KatgertI (|2006l ) for a subsample of 54 galaxies 
finding complete agreement. 

Adopting a weighted fit (see Eq.[8l we derived the FP: 
re oca-^ °'*(7)J°-*^ for the 141 early-type R<18 galaxies. The 
low value of a can be related to the extension of our sam- 
ple towards very low values of velocity dispersion and and 
also to the method used to determine a. Considering the 
subsample of ctq > 100 km s~^ galaxies the FP turns out 
to follow the relation rcOa a^''^^ (I)'^'^'^^ which is consistent 
with the Coma FP obtained by JFK96. Observing the sig- 
nificant change in the a value between the total and the 
high-ao samples, we investigate the possible curvature of 
the FP. Applying the same cuts to a simulated sample fol- 
lowing a linear FP relation, we ascertain that the observed 
curvature can be explained by selection effects. Departure 
from a linear trend of the FP relation has been claimed by 
D'Onofrio et al. (2008), and Desroches et al. (2001) found 
that the faint-end luminosity cut influences the coefficient a 
of the FP, in agreement with Nigoche-Netro (2008). Similar 
results were found by Hyde & Bernardi (2008) studying the 
FP residuals along the plane. However, these works pointed 
out that conclusive evidence of the FP curvature needs ei- 
ther more robust statistics or higher accuracy in the velocity 
dispersion estimates for low mass systems. 

The most important result of this work is our demon- 
stration that the FP residuals are correlated with stellar 
population characteristics (line- strength indices and de- 
rived SSP parameters). In particular, FP residuals are anti- 
correlated both with the a-element abundance ratio, a/Fe, 
and with galaxy age resulting in trends whereby galaxies 
more compact than expected from the FP relation have stel- 
lar pupulations systematically older and with higher abun- 
dances than average. Previous studies have reported corre- 
lations of the FP residuals with stellar age (e.g. Forbes et 
al. 1998; Reda et al. 2005). Ahhough our FP residual do 
show a correlation with age, we recover a much stronger sig- 
nal for a/Fe than for age, suggesting that this is the more 
fundamental dependence. Indeed, a multiple regression anal- 
ysis suggests there is no age correlation at fixed a/Fe, but 
a very strong a/Fe correlation at fixed age. The correlation 
between a/Fe and FP residuals does not indicate a direct 
causal effect, since varying a/Fe at fixed age and Z/H has 
little effect on the stellar mass-to-light. Instead, the corre- 
lation suggests that the structural properties and the star- 
formation history are both dependent on some unobserved 
aspect of the galaxy assembly process. At face value, our re- 
sults are broadly consistent with recent galaxy merger sim- 
ulations, which predict a sequence of formation mechanisms 
governed by the varying importance of dissipation (Hopkins 
et al. 2008). In this scenario, mergers with a higher initial 
gas-fraction trigger more centrally-concentrated starbursts, 
and higher a abundances due to the short duration of star- 
formation in the burst. 



ACKNOWLEDGMENTS 

AG gratefully acknowledges the hospitality of the Univer- 
sity of Birmingham during her stays there, where some of 
the work was performed. CPH acknowledges financial sup- 
port from STFC. RJS is supported by STFC rolling grant 
PP/C501568/1 "Extragalactic Astronomy and Cosmology 
at Durham 2005-2010". 

This work was carried out in the framework of the col- 
laboration of to the FP7-PEOPLE-IRSES-2008 project AC- 
CESS "A Complete CEnsus of Star-formation and nuclear 
activity in the Shapley supercluster" . 

We thank the anonymous referee for his/her comments 
which helped to improve this work. 



REFERENCES 

Beers T. C, Flynn K., Gebhardt K., 1990, AJ, 100, 32 
Bender R., Burstein D., Faber S. M., 1992, ApJ, 399, 462 
Bernardi M., et al., 2003, AJ, 125, 1866 
Berlin E., Arnouts S., 1996, A&A, 117, 393 
Berlin G., Ciotti L., Del Principe M., 2002, A&A, 386, 149 
Binney & Merrifield, 1998, in Galactic Astronomy, Prince- 
ton Univeristy Press 
Caon, N., Capaccioli M., D'Onofrio M., 1993, MNRAS, 
265, 1013 

Cappellari M., et al., 2006, MNRAS, 366, 1126 
Djorgovski S., Davis M., 1987, ApJ, 313, 59 
D'Onofrio et al, 2008, arXiv:0804. 18921 
Dressier et al., 1987, ApJ, 313, 42 
Faber S. M. & Jackson R. E., 1976, ApJ, 204, 668 
Forbes et al., 1998, ApJ, 508, 43 

Fukugita M., Shimasaku K., Ichikawa T. A. W., 1995, 

PASP, 107, 945F 
Gerhard O., Kronawitter A., Sagha R. P., Bender R. 2001, 

AJ, 121, 1936 
Graham A. W., Guzman R., 2003, AJ, 125, 2936 
Gregg M., 1992, ApJ, 384, 43 
Greggio L & Renzini A., 1983, A. A, 118, 217. 
Guzman R. & Lucey J. R., 1993, MNRAS, 263, L47. 
Guzman R., Lucey J. R., Bower R. G. 1993, MNRAS, 265, 

731. 

Haines C. P., Merluzzi P., Mercurio A., Gargiulo A., Kru- 

sanova N., Busarello G., La Barbera F., Capaccioli M., 

2006, MNRAS, 371, 55 
Haines C. P., Gargiulo A., La Barbera F., Mercurio A., 

Merluzzi P., Busarello G., 2007, MNRAS, 381, 7 
Hopkins P. F., Cox T. J., Hernquist L., 2008, ApJ, 689, 17 
j0rgensen I., Franx M., 1994, ApJ, 433, 553 
j0rgensen I., Franx M., Kjaegaard P., 1995, MNRAS, 276, 

1341 

J0rgensen I., Franx M., Kjaegaard P., 1996, MNRAS, 280, 
167 (JFK96) 

Kelson D. D., lUingworth G. D., van Dokkum P. G., Franx 

M., 2000, ApJ, 531, 184 
Kormendy J., 1976, ApJ, 218, 333 
Kroupa P., 2001, MNRAS, 322, 231 
Kuntschner et al., 2006, MNRAS, 369, 497 
La Barbera F., Busarello G., Merluzzi P., de la Rosa I., 

Coppola G., Haines C. P., 2008a, ApJ, 689, 913 
La Barbera F., de Carvalho R. R., Kohl-Moreira J. L., 



16 Gargiulo et al. 



Gal R. R., Soares-Santos M., Capaccioli M., Santos R., 

Sant'Anna N., 2008b, PASP, 120, 681 
Maraston C, 2005, MNRAS, 362, 799 
Matkovic A. Guzman R., 2005, MNRAS, 362, 289 
Matteucci F., 1996, Fundamentals of Cosmic Physics, Vol. 

17, pp. 283-396 
Mercurio A., Merluzzi P., Haines CP., Gargiulo A., Kru- 

sanova N., Busarello G., La Barbera F., Capaccioli M., 

Govone G., 2006, MNRAS, 368, 109 
Nelan J. E., Smith R. J., Hudson M. J., Wegner G. A., 

Lucey J. R., Moore S. A. W., Quinney S. J., Suntzeff N. 

B., 2005, ApJ, 632, 137 
Nigoche-Netro A. et al., 2008, preprint, asrto-ph 08051142 
Pahre M. A., Djorgovski S. G., de Garvalho R. R., 1999, 

ASP Conference Series, Vol. 163, 17 
Poggianti B. M., 1997, A&AS, 122, 399 
Proctor R. N. & Sansom A. E., 2002, MNRAS, 333, 517 
Prugniel P. & Simien F., 2002, A&A, 309, 749 
Rawle T. D., Smith R. J., Lucey J. R., Swinbank A., 2008, 

MNRAS, 389, 1891 
Reda F. et al., 2005, MNRAS, 360, 693 
Robertson B., Cox T. J., Henquist L., Franx M., Hopkins 

P. F., Martini P., Springel V., 2006, ApJ, 641, 21 
Sanchez-Blazquez P., et al., 2006, MNRAS, 371, 703 
Sanchez-Blazquez P., Forbes D. A.; Strader J., Brodie J., 

Proctor R., 2007, MNRAS, 377, 759 
Schlegel D. J. et al., 1998, ApJ, 500, 525 
Smith R. J., et al., 2004, AJ, 128, 1558 
Smith R. J., Lucey J. R., Hudson M. J., 2007, MNRAS, 

381, 1035 (SLH) 
Smith R. J., Lucey J. R., Hudson M. J. 2009, MNRAS, 

submitted 

Steinmetz M. & Navarro J. F., 2002, NewA, 7, 155 
Thomas D., Maraston C, Bender R., 2003, MNRAS, 339, 
897 

Thomas D., Maraston C, Korn A., 2004, MNRAS, 351, 
L19 

Thomas J., Saglia R. P., Bender R., Thomas D., Gebhardt 
K., Magorrian J., Corsini E. M., Wegner G., 2005, MN- 
RAS, 360, 1355 

Thomas T., Katgert P., 2006, A&A, 446, 31 

Tinsley B. M., 1980, Fund Cosmic Phys., 5, 287 

Toomre A., Toomre J., 1972, ApJ, 178, 623 

Toomre A., 1977, in Evolution of Galaxies and Stellar Pop- 
ulations, ed. B. M. Tinsley & R. B. Larson, 401 

Trujillo L, Burkert A., BeU E. F., 2004, ApJL, 600, 39 

Vandame B., 2004, PhD thesis 

Vazdekis A., Trujillo L, Yamada Y. 2004, ApJ, 601L 
Young C. K., Currie M. J., 1994, MNRAS, 268, Lll 
Wainer H. & Thissen D., 1976, Psychometrika, 41, 9 
Wuyts S. et al., 2004, ApJ, 605, 677 

Zaritsky D., Gonzalez A. H., Zabludoff A. L, 2006, ApJ, 
638, 725 



